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Abstract 


Recent  work  at  the  Defence  Research  Establishment  Atlantic  (DREA)  has  led  to 
the  development  of  software  for  predicting  the  target  strength  of  underwater  ob¬ 
jects.  This  capability  is  incorporated  into  the  general  acoustic  radiation  predic¬ 
tion  software  tool,  AVAST.  While  AVAST  has  been  successfully  validated  against 
a  limited  set  of  theoretical  data,  it  has  not  been  fully  validated  against  experimental 
data  or  other  numerical  codes.  In  this  report,  AVAST  predictions  are  compared  to 
those  made  by  an  axisymmetric-based  code  (ASBIEM)  developed  at  SACLANT- 
CEN  for  small  cylinders  with  either  flat  or  hemispherical  endcaps.  The  results  were 
compared  at  three  frequencies  for  the  full  monostatic  case  and  the  bistatic  case 
(with  three  different  source  locations).  For  the  bistatic  case,  there  was  generally 
good  agreement  between  AVAST,  ASBIEM,  and  theoretical  results  for  the  specific 
source  locations  selected.  The  results  from  the  cylinder  with  hemispherical  end- 
caps  showed  slightly  better  agreement  than  those  from  the  flat  end  cylinder.  For 
the  monostatic  analyses,  the  AVAST  and  ASBIEM  codes  showed  very  good  agree¬ 
ment  at  the  lower  frequencies  for  both  cylinder  models.  At  the  lowest  frequency, 
AVAST  seemed  to  show  a  slightly  better  agreement  with  the  theoretical  results.  At 
the  highest  frequency,  however,  AVAST,  while  correctly  matching  the  broadside 
and  endfire  results  and  the  general  level,  appears  to  contain  some  anomalous  side 
lobes  and  nulls.  In  general,  the  comparisons  showed  that  AVAST  and  ASBIEM 
gave  comparable  results  when  used  to  predict  rigid  target  strength.  The  ASBIEM 
code  runs  on  the  order  of  one  magnitude  faster  than  the  AVAST  code,  but  cannot  be 
used  to  model  structures  which  are  not  axisymmetric. 


Resume 


Des  travaux  rdcents  effectues  au  Centre  de  recherches  pour  la  Defense  Atlantique 
(CRD A)  ont  conduit  k  1’ Elaboration  d’un  logiciel  pour  la  prevision  de  l’intensitd 
de  cibles  sous-marines.  Cette  capacity  est  incorporde  k  l’outil  logiciel  general  de 
provision  du  rayonnement  acoustique,  AVAST.  Ce  logiciel  a  dtd  valide  avec  succds 
pour  un  ensemble  limitd  de  donndes  thdoriques,  mais  il  n’a  pas  dtd  entierement 
validd  pour  des  donndes  experimentales  ou  pour  d’autres  codes  numdriques.  Dans 
ce  rapport,  les  prdvisions  fournies  par  AVAST  sont  compardes  aux  prdvisions  obtenues 
avec  un  logiciel  basd  sur  des  structures  axisymdtriques  (ASBIEM)  dlabore  au  SACLANT- 
CEN  pour  de  petits  cylindres  k  bouchons  d’extrdmitds  plats  ou  hdmisphdriques.  On 
a  compard  les  rdsultats  obtenus  k  trois  frdquences  pour  le  cas  entidrement  monos- 
tatique  et  pour  le  cas  bistatique  (avec  trois  emplacements  diffdrents  de  la  source). 

Pour  le  cas  bistatique,  1’ accord  etait  gdndralement  bon  entre  les  rdsultats  donnes 
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par  AVAST,  les  r6sultats  donnas  par  ASBIEM  et  les  resultats  theoriques  pour  les 
emplacements  choisis  de  la  source.  L’  accord  des  resultats  obtenus  avec  le  cylindre 
k  bouchons  d’extremites  h£mispheriques  etait  16g£rement  superieur  a  celui  obtenu 
avec  le  cylindre  a  bouchons  d’ extremity  plats.  Dans  le  cas  des  analyses  monosta- 
tiques,  les  logiciels  AVAST  et  ASBIEM  ont  donne  un  tr&s  bon  accord  aux  frequences 
inferieures  pour  les  deux  types  de  cylindres.  A  la  plus  basse  frequence,  AVAST 
semblait  donner  un  accord  l<§gferement  superieur  avec  les  resultats  theoriques.  A  la 
plus  haute  frequence,  cependant,  AVAST,  raeme  s’il  donnait  des  resultats  en  accord 
avec  ceux  obtenus  pour  un  rayonnement  transversal  et  un  rayonnement  longitudi¬ 
nal,  de  meme  que  pour  le  niveau  general,  semblait  contenir  certains  lobes  lateraux 
et  creux  anormaux.  En  general,  les  comparaisons  ont  montre  que  les  logiciels 
AVAST  et  ASBIEM  donnaient  des  resultats  comparables  lorsqu’ils  etaient  utilises 
pour prevoir  l’intensite des  cibles  rigides.  L’ execution  du  logiciel  ASBIEM est  d’un 
ordre  de  grandeur  plus  rapide  que  celle  du  logiciel  AVAST,  mais  elle  ne  permet  pas 
de  modeliser  les  structures  non  axisymetriques. 


ii 
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Executive  summary 


Introduction 

An  important  piece  of  information  in  the  detection  and  localization  of  underwater 
objects  is  their  acoustic  target  strength.  The  likely  targets  can  range  in  size  from 
ships  and  submarines  to  sea  mines  over  a  broad  range  of  frequencies.  For  sim¬ 
ple  geometric  shapes  (e.g.,  spheres,  infinite  cylinders)  and  structural  compositions 
(e.g.,  rigid  walls)  theoretical  target  strength  solutions  exist  and  are  used  to  approx¬ 
imate  the  behaviour  of  real  targets.  However,  these  approximations  are  of  limited 
use  if  the  object  is  sufficiently  complex.  For  some  targets,  such  as  submarines, 
some  experimental  data  exist,  but  the  measurements  are  rarely  applicable  across 
all  submarines.  Thus,  numerical  methods  which  can  be  used  over  a  wide  range 
of  structures,  frequencies,  and  environmental  conditions  to  predict  acoustic  target 
strength  are  valuable  tools  for  the  sonar  analyst.  A  variey  of  numerical  methods  are 
available  with  most  being  based  on  either  the  finite  element  or  boundary  element 
method. 

Recent  work  at  DREA  has  led  to  the  development  of  a  mixed  finite  element  and 
boundary  element  acoustic  radiation  prediction  software  tool,  AVAST,  which  can 
be  used  to  predict  the  acoustic  target  strength  of  full  three-dimensional  underwa¬ 
ter  structures.  Also,  in  a  separate  program  at  SACLANTCEN,  an  axisymmetric 
boundary  element-based  code,  ASBIEM,  was  developed  to  investigate  the  rigid  tar¬ 
get  strength  characteristics  of  mine-like  targets  over  a  large  frequency  range.  In  the 
work  done  at  SACLANTCEN,  rigid  target  strength  predictions  were  made  for  two 
mine-like  targets,  each  a  cylinder  of  overall  length  of  2.0m  and  radius  of  0.25m, 
one  with  flat  endcaps  and  one  with  hemispherical  endcaps.  In  this  technical  mem¬ 
orandum,  the  theoretical  background  to  the  two  computer  codes  are  discussed  and 
then  both  methods  are  used  to  predict  the  target  strength  for  the  small  2m  cylinders 
at  three  different  frequencies  with  a  variety  of  source  locations  and  the  results  are 
compared  to  each  other  and  to  selected  theoretical  results.  Results  were  generated 
for  both  the  bistatic  and  full  monostatic  cases. 


Principal  Results 

For  the  bistatic  case,  there  was  generally  good  agreement  between  AVAST,  AS¬ 
BIEM,  and  theoretical  results  for  the  specific  source  locations  selected.  The  results 
from  the  cylinder  with  hemispherical  endcaps  seem  to  show  better  agreement  than 
those  from  the  flat  end  cylinder.  The  dominant  feature  of  the  bistatic  target  strength 
in  general  was  forward  scattering  which  was  always  the  largest  amplitude  lobe. 

For  the  monostatic  analyses,  the  AVAST  and  ASBIEM  codes  showed  very  good 
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agreement  for  both  the  2.5  kHz  and  5  kHz  cases  for  both  cylinder  models.  At  the 
lowest  frequency,  AVAST  seemed  to  show  a  slightly  better  agreement  with  the  the¬ 
oretical  results.  At  the  highest  frequency,  however,  AVAST,  while  correctly  match¬ 
ing  the  broadside  and  endfire  results  and  the  general  level,  appears  to  contain  some 
anomalous  side  lobes  and  nulls. 


Significance  of  Results 

In  general,  the  above  comparisons  show  that  AVAST  and  ASBIEM  give  comparable 
results  when  used  to  predict  rigid  target  strength.  The  ASBIEM  code  runs  on  the 
order  of  one  magnitude  faster  than  the  AVAST  code,  but  cannot  be  used  to  model 
structures  which  are  not  axisymmetric.  At  the  highest  frequency  in  the  monostatic 
case,  AVAST  predicted  a  few  apparently  anomalous  side  lobes  and  nulls.  In  past 
analyses,  these  side  lobes  disappear  when  a  more  refined  grid  is  used,  so  while  the 
bistatic  results  implied  the  grid  was  sufficiently  refined,  these  results  indicate  this 
was  not  the  case.  At  least  for  the  AVAST  code,  ten  elements  per  acoustic  wavelength 
is  not  quite  sufficient  for  the  complete  analysis. 


Future  Plans 

AVAST  supports  a  number  of  capabilities  not  explored  in  this  report  which  may 
be  of  further  interest  to  those  involved  in  submarine  detection  or  mine  counter¬ 
measures.  Issues  such  as  target  strength  prediction  in  shallow  water  including  the 
effects  of  the  bottom  and  the  water  surface,  and  target  strength  of  partially  buried 
mines  can  be  examined  using  the  existing  software.  Future  work  on  AVAST  will 
concentrate  on  such  things  as  examining  elastic  target  strength  (allowing  for  inter¬ 
nal  structure),  exploring  ways  to  increase  the  frequency  range  (including  examining 
the  anomalous  lobes  demonstrated  here),  and  examining  a  variety  of  environmental 
conditions  such  as  shallow  water  and  varying  bottom  conditions. 


L.E.  Gilroy,  J.  Fawcett;  2001 ;  Comparison  of  Three-Dimensional  and 
Axisymmetric  Software  for  Predicting  Acoustic  Target  Strength; 
DREATM  2001-018;  Defence  Research  Establishment  Atlantic. 
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Sommaire 


Introduction 

Un  element  d’information  important  dans  la  detection  et  la  localisation  de  cibles 
sous-marines  est  leur  intensity  acoustique.  Les  cibles  probables  peuvent  avoir  des 
tailles  allant  de  celle  des  navires  et  des  sous-marins  k  celle  des  mines  marines  sur 
une  vaste  gamme  de  frequences.  Pour  les  formes  geometriques  simples  (p.  ex. 
spheres,  cylindres  infinis)  et  les  elements  de  structure  (p.  ex.  murs  rigides)  il  existe 
des  solutions  pour  l’intensite  des  cibles  theoriques  et  elles  sont  utilisees  pour  obtenir 
une  approximation  du  comportement  des  cibles  rdelles.  Cependant,  ces  approxima¬ 
tions  sont  d’ utilisation  limitee  lorsque  l’objet  est  assez  complexe.  Pour  certaines 
cibles,  par  exemple  les  sous-marins,  il  existe  certaines  donndes  experimentales, 
mais  les  mesures  sont  rarement  applicables  a  tous  les  sous-marins.  Par  consequent, 
les  methodes  numeriques  qui  peuvent  etre  utilisees  pour  une  vaste  gamme  de  struc¬ 
tures,  frequences  et  conditions  environnementales  en  vue  de  prevoir  l’intensite 
acoustique  des  cibles  sont  des  outils  prdcieux  pour  le  specialiste  de  l’analyse  sonar. 
Une  varidte  de  methodes  numeriques  sont  disponibles,  et  la  plupart  sont  basees  soit 
sur  la  methode  par  elements  finis,  soit  sur  la  methode  par  elements  frontieres. 

Des  travaux  recents  du  CRDA  ont  conduit  k  l’eiaboration  d’un  outil  logiciel  de 
prevision  du  rayonnement  acoustique  combinant  la  methode  par  elements  finis  et  la 
methode  par  elements  frontieres,  AVAST,  qui  permet  de  prevoir  l’intensite  acous¬ 
tique  de  structures  sous-marines  entierement  tridimensionnelles.  De  plus,  dans 
le  cadre  d’un  programme  distinct  mend  au  SACLANTCEN,  un  logiciel  base  sur 
des  elements  frontieres  axisymdtriques,  ASBIEM,  a  dtd  dlabord  pour  1’dtude  des 
caractdristiques  d’intensitd  de  cibles  rigides  semblables  a  des  mines  sur  une  vaste 
gamme  de  frequences. 

Dans  le  cadre  des  travaux  effectuds  au  SACLANTCEN,  des  previsions  de  l’intensitd 
de  cibles  rigides  ont  dtd  faites  pour  deux  cibles  semblables  k  des  mines,  consistant 
chacune  en  un  cylindre  ayant  une  longueur  hors  tout  de  2,0  m  et  un  rayon  de  0,25  m, 
une  k  bouchons  d’extrdmitds  plats,  1’ autre  a  bouchons  d’extrdmitds  hdmisphdriques. 
Dans  ce  document  technique,  on  traite  du  contexte  thdorique  des  deux  logiciels,  puis 
on  utilise  les  deux  mdthodes  pour  prdvoir  l’intensitd  des  petits  cylindres  cibles  de  2 
m  k  trois  frequences  diffdrentes  avec  une  varidte  d’ emplacements  de  la  source  et  on 
compare  les  rdsultats  entre  eux  et  k  des  rdsultats  thdoriques  choisis.  Des  rdsultats 
ont  dtd  obtenus  pour  le  cas  bistatique  et  le  cas  entidrement  monostatique. 
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Principaux  resultats 

Pour  le  cas  bistatique,  V accord  etait  gen^ralement  bon  entre  les  resultats  donnas 
par  AVAST,  les  resultats  donnes  par  ASBIEM  et  les  resultats  theoriques  pour  les 
emplacements  choisis  de  la  source.  L’  accord  des  resultats  obtenus  avec  le  cylindre 
a  bouchons  d’extr^mitds  h6misph6riques  etait  l^g^rement  superieur  a  celui  obtenu 
avec  le  cylindre  k  bouchons  d’extremit6s  plats.  La  caracteristique  dominante  de 
l’intensite  de  cible  bistatique  en  general  etait  la  diffusion  vers  l’avant,  qui  etait 
toujours  le  lobe  de  plus  grande  amplitude. 

Dans  le  cas  des  analyses  monostatiques,  AVAST  et  ASBIEM  ont  donne  un  tr&s 
bon  accord  k  2,5  kHz  et  k  5  kHz  pour  les  deux  types  de  cylindres.  A  la  plus 
basse  frequence,  AVAST  semblait  donner  un  accord  legSrement  superieur  avec  les 
resultats  theoriques.  A  la  plus  haute  frequence,  cependant,  AVAST,  raeme  s’il  don- 
nait  des  resultats  en  accord  avec  ceux  obtenus  pour  un  rayonnement  transversal  et 
un  rayonnement  longitudinal,  de  meme  que  pour  le  niveau  general,  semblait  con- 
tenir  certains  lobes  lateraux  et  creux  anormaux. 

Signification  des  resultats 

En  general,  les  comparaisons  susmentionnees  montrent  que  les  logiciels  AVAST 
et  ASBIEM  donnent  des  resultats  comparables  lorsqu’ils  sont  utilises  pour  prevoir 
l’intensite  de  cibles  rigides.  L’ execution  du  logiciel  ASBIEM  est  d’un  ordre  de 
grandeur  plus  rapide  que  celle  du  logiciel  AVAST,  mais  elle  ne  permet  pas  de 
modeiiser  les  structures  non  axisymetriques.  A  la  plus  haute  frequence  du  cas 
monostatique,  AVAST  prdvoyait  quelques  lobes  et  creux  apparemment  anormaux. 
Dans  les  analyses  anterieures,  ces  lobes  lateraux  disparaissaient  lorsqu’on  utilisait 
un  quadrillage  plus  perfectionne.  Par  consequent,  meme  si  les  resultats  du  cas  bis¬ 
tatique  laissaient  sous-entendre  que  le  quadrillage  etait  suffisamment  perfectionne, 
les  nouvelles  donnees  indiquent  qu’il  ne  l’etait  pas.  Pour  AVAST  au  moins,  il  n’est 
pas  tout  a  fait  suffisant  d’utiliser  dix  elements  par  longueur  d’onde  acoustique  pour 
1’  analyse  complete. 


Plans  futurs 

AVAST  permet  un  certain  nombre  de  fonctions  qui  ne  sont  pas  etudiees  dans  ce 
rapport  et  qui  peuvent  interesser  davantage  les  responsables  de  la  detection  des 
sous-marins  et  de  la  lutte  contre  les  mines.  Les  questions  telles  que  la  prevision 
de  l’intensite  des  cibles  en  eaux  peu  profondes,  y  compris  les  effets  du  fond  et  de 
la  surface  de  l’eau,  et  de  l’intensite  des  mines  partiellement  enfouies  peuvent  etre 
etudiees  a  1’aide  du  logiciel  actuel.  Les  travaux  futurs  relatifs  a  AVAST  seront  axes 
sur  des  aspects  tels  que  l’etude  de  l’intensite  des  cibles  eiastiques  (tenant  compte  de 


vi 


DREA  TM  2001-018 


P516761.PDF  [Page:  12  of  37] 


la  structure  interne),  la  recherche  de  moyens  pour  accroitre  la  plage  de  frequences  (y 
compris  l’etude  des  lobes  anormaux  observes)  et  l’etude  d’une  variety  de  conditions 
environnementales,  par  exemple  les  conditions  d’eaux  peu  profondes  et  de  fonds 
variables. 


L.E.  Gilroy,  J.  Fawcett;  2001 ;  Comparison  of  Three-Dimensional  and 
Axisymmetric  Software  for  Predicting  Acoustic  Target  Strength; 
DREA  TM  2001-018;  Centre  de  recherches  pour  la  defense 
Atlantique. 
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1  Introduction 


The  ability  to  predict  the  target  strength  of  underwater  objects  is  of  great  interest  to 
the  sonar  community  to  assist  in  target  localization  and  identification.  Numerical 
methods  of  target  strength  prediction  are  useful  as  experimental  measurements  of 
target  strength  are  difficult  to  perform  and  cannot  easily  cover  the  range  of  condi¬ 
tions  available  to  numerical  methods.  A  variey  of  numerical  methods  are  available 
with  most  being  based  on  either  the  finite  element  or  boundary  element  method. 

Recent  work  in  support  of  the  Ship  Noise  Project  at  the  Defence  Research  Estab¬ 
lishment  Atlantic  (DREA)  has  led  to  the  development  of  a  mixed  finite  element  and 
boundary  element  acoustic  radiation  prediction  software  tool,  AVAST  [1,  2,  3,  4], 
This  software  originally  was  used  to  predict  the  low  frequency  radiated  noise  from 
submerged  or  floating  elastic  structures  and  is  intended  to  assist  DND  in  either 
optimizing  the  structural  arrangement  to  minimize  radiated  noise  or  in  examining 
existing  structures  to  isolate  noise  sources.  During  the  development  of  the  soft¬ 
ware,  a  capability  to  predict  scattered  sound  pressure  (given  an  external  source) 
was  also  implemented.  This  capability  allowed  for  both  a  rigid  structure  and  an 
elastic  structure  assuming  that  a  natural  frequency  analysis  of  the  structure  could 
be  performed  up  to  the  required  frequency.  Given  the  existing  scattering  capabil¬ 
ity,  a  target  strength  prediction  capability  was  a  relatively  simple  addition  to  the 
AVAST  software.  As  with  the  scattering  problem,  AVAST  can  solve  both  the  rigid 
and  elastic  target  cases,  however  further  investigation  is  required  to  determine  if  a 
sufficiently  large  range  of  frequencies  can  be  accommodated  for  the  elastic  target 
case.  AVAST  has  a  full  3D  capability,  so  can  it  be  used  to  examine  complex  asym¬ 
metric  structures.  However,  AVAST  does  not  have  an  axisymmetric  capability  and, 
thus,  requires  a  much  larger  model  for  such  targets  than  would  be  required  with  a 
purely  axisymmetric  code. 

In  support  of  the  project  on  Modelling  Scattering  from  Buried  Objects,  an  axisym¬ 
metric  boundary  element-based  (ASBIEM)  program,  was  developed  at  SACLANT- 
CEN  [5].  The  program  was  used  to  investigate  the  rigid  and  elastic  target  strength 
characteristics  of  mine-like  targets  over  a  large  frequency  range.  The  advantage  of 
such  a  program  is  the  capability  to  cover  a  very  large  frequency  range  with  rela¬ 
tively  few  elements  (as  one  is  only  required  to  model  the  2D  shape  of  the  target 
using  line  elements).  The  primary  disadvantage  of  such  a  program  is  the  restriction 
to  axisymmetric  geometries. 

Initial  target  strength  work  with  AVAST  has  involved  the  prediction  of  rigid  tar¬ 
get  strength  for  a  variety  of  targets  focussed  on  submarine-based  models,  including 
large  cylinders  and  more  realistic  submarine  shapes  [6,  7,  8].  In  the  work  done 
at  SACLANTCEN,  rigid  target  strength  predictions  were  made  for  two  mine-like 
targets,  each  a  cylinder  of  overall  length  of  2.0m  and  radius  of  0.25m,  one  with 
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flat  endcaps  and  one  with  hemispherical  endcaps.  In  this  technical  memorandum, 
the  theoretical  background  to  the  two  computer  codes  are  discussed  and  then  both 
methods  are  used  to  predict  the  target  strength  for  the  small  2m  cylinders  at  three 
different  frequencies  (2.5  kHz,  5.0  kHz,  and  7.5  kHz)  with  a  variety  of  source  loca¬ 
tions  and  the  results  are  compared  to  each  other  and  to  selected  theoretical  results. 


2  AVAST 


Under  contract  to  DREA,  Martec  Ltd.  has  developed  a  series  of  computer  programs, 
collectively  named  AVAST,  for  use  in  the  numerical  prediction  of  the  acoustic  radia¬ 
tion  and  scattering  from  floating  or  submerged  elastic  structures  immersed  in  either 
infinite,  half-space  or  finite  depth  fluid  domains.  More  recently,  the  modelling  ca¬ 
pabilities  of  the  AVAST  code  have  been  extended  to  include  target  strength  analysis 
of  both  elastic  and  inelastic  structures.  AVAST  combines  both  the  finite  element 
method  for  the  structure  and  the  boundary  integral  equation  technique  for  the  fluid. 
Note  that,  for  a  rigid  target,  a  structural  finite  element  model  is  not  required,  it  is 
simply  necessary  to  model  the  shape  of  the  wetted  surface  of  the  structure. 


For  problems  involving  exterior  fluid  domains,  AVAST  employs  a  boundary  ele¬ 
ment  formulation  based  upon  the  classical  Helmholtz  integral  equation,  i.e.: 


a 


c,(Pf)p(P/)=  /  G(P,,  Q)|Pa,!»(0) 


9G(P/,Q) 

dno 


p(Q ) 


dS{Q)  +  Airpt(Pf) 


(1) 

where  G  represents  the  Green’s  function  due  to  a  harmonic  point  source  located  at 
point  Q,  p{Pf)  represents  the  acoustic  pressure  at  a  point  Pf  in  the  fluid  domain  (sat¬ 
isfying  the  Helmholtz  differential  equation  for  time-harmonic  waves,  V2p(Pf)  + 
k2p(Pf)  =  0),  u>  represents  the  angular  frequency,  p  represents  the  fluid  density, 
u(Q)  represents  the  displacement  amplitude  normal  to  the  body  surface  S  at  the 
point  Q,  Pi{Pf)  represents  the  prescribed  pressure  of  the  incident  acoustic  wave, 
and  7 iq  represents  the  surface  normal  vector  at  point  Q. 


For  the  special  case  of  a  body  submerged  in  an  infinite  fluid  domain,  this  Green’s 
function  can  be  expressed  using  the  following  formula: 

-ikR(P,,Q) 

=  (2) 

where  R  is  the  distance  between  Pf  and  Q  and  k  represents  the  acoustic  wavenum¬ 
ber  (t o/c),  where  c  represents  the  fluid  sound  speed.  In  addition,  the  coefficient 
Cg(Pf)  is  a  function  of  the  local  geometry,  having  a  value  of  47r  for  locations  Pf 
fully  immersed  in  the  fluid  domain  and  a  value  of  2 ir  for  points  residing  on  the 
body  surface. 
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In  order  to  solve  the  integral  formulation  posed  above  in  Equation  1,  a  four-step  nu¬ 
merical  technique  was  developed  and  implemented  into  the  AVAST  code.  The  first 
step  involves  the  discretization  of  the  wet  surface  of  the  underwater  target  into  a  col¬ 
lection  of  either  3-node  or  4-node  boundary  element  panels.  Step  two  involves  the 
numerical  integration  of  Helmholtz  integral  equation  over  the  surface  of  the  body 
(represented  by  the  boundary  element  panels).  This  leads  to  a  system  of  algebraic 
equations  relating  the  acoustic  pressure  on  the  surface  of  the  body  to  the  incident 
pressure  field  and  surface  normal  velocities,  i.e.: 


[H]  {p}  =  [G]  {«}  +  {pj 


(3) 


where  the  number  of  unknown  surface  pressures  (length  of  the  vector  {p})  is  equal 
to  the  number  of  boundary  element  panels  and  the  matrix  operators  [H]  and  [G]  can 
be  expressed  using  the  following  forms: 


and 


(4) 

(5) 


where  the  subscripts  ij  represent  the  z-th  row  and  the  j-th  column  of  the  matrices 
[H]  and  [G], 

The  third  step  in  the  AVAST  solution  procedure  involves  solving  Equation  3  for  the 
surface  acoustic  pressures,  i.e.: 


{p}  =  [ff]-1(^3[G]W  +  fe})  (6) 


These  surface  acoustic  pressures  can  also  be  expressed  in  terms  of  nodal  forces, 
{/ a }  (on  the  structure),  i.e.: 

if  a}  =  [T]T{A}{Pi} 

=  pu>>  [If  [A]  [H]-1  [G]  [T]  {x}  +  [If  [A]  [H]~l  {Pi} 

=  [Mf]  {x}  +  [if  [A]  [. H [-1  {pj  (7) 


where  [T]  represents  a  transformation  matrix  relating  the  global  cartesian  structural 
degrees  of  freedom  to  the  boundary  element  panel  normals,  [A]  represents  a  diago¬ 
nal  fluid  panel  area  matrix,  [Mf]  represents  the  so-called  acoustic  fluid  mass  matrix, 
and  {x}  represents  a  vector  of  body  displacements  defined  in  terms  of  the  global 
coordinate  system. 


In  the  fourth  and  final  stage  of  an  AVAST  solution  technique  the  surface  acoustic 
pressures,  computed  during  the  third  phase,  are  used  to  predict  the  sound  pressure 
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levels  at  points  residing  in  the  fluid  domain.  Once  these  field  pressures  have  been 
computed,  estimates  of  the  body’s  target  strength  can  be  generated  using  Equa¬ 
tion  16. 

AVAST  models  can  be  imported  from  a  number  of  sources  including  some  standard 
mesh  generation  packages.  Using  the  AVAST  GUI,  sources,  field  points,  wave 
guides,  free  surfaces,  fluid  properties,  etc.  may  be  added  to  the  model  which  is  then 
solved  at  specified  frequencies  or  ranges  of  frequencies.  Results  may  be  plotted  in 
the  AVAST  GUI  or  exported. 


3  ASBIEM 


The  azimuthally-symmetric  Boundary  Integral  Equation  Method  (ASBIEM)  starts 
from  the  three-dimensional  integral  equation  for  the  pressure  on  a  rigid  object,  fol¬ 
lowing  the  form  of  Eq.(l)  (with  velocities,  u,  set  to  zero  as  it  is  a  rigid  object) 


1  fdG(rf,Zf,rQ,zQ,ef-dQ) 

P  =  ~^JS - dn^ - dS(Q)  +  p. 


(8) 


Here  a  cylindrical  coordinate  system  is  used;  the  object  is  assumed  to  be  azimuth- 
ally-symmetric  about  the  z-axis;  the  angle  6  is  the  polar  angle  around  the  z-axis. 
The  Green’s  function  itself  depends  only  on  the  difference  of  the  field  angle  0f 
and  the  integration  angle  Oq.  To  exploit  the  angular  symmetry  of  the  problem,  the 
pressure  field  p(r/5  z/,  9f)  is  expanded  in  a  Fourier  series  with  unknown  coefficients 
Pm(r,z), 

OO 

p{rf,zf,0f)  =  ^2  pm(rf,Zf)e,Tne.  (9) 

Tfl  —  — OO 

Similarly, 

OO 

P*=  ^2  i”*(r7>*/)eimS  (1°) 


where  in  this  case  the  incident  coefficients  are  known.  The  Green’s  function  or  its 
normal  derivative  can  also  be  expanded  in  the  form. 


where 


8G(P,,Q ) 

dnq 


OO 

Y.  9m{rf,Zf,rQ,zq)eimT 

771—  —  OO 


9m{rf,zf,rq, 


j  zfiTQi  zqiT )eimrdT 


(ID 


(12) 
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andr  =  9f-6Q.  In  fact  since  G(rf,zf,rQ,zQ,T)  =  G(rf,zf,rQ,zQ,~T),  Eq(12) 
can  be  rewritten  as 

9m{rf,zf,rq,zq)  =  2  /  G{rf,  zf,  rQ,  zq,  r)  cos(mr)(ir.  (13) 

Jo 

Using  the  series  expressions  in  Eq.(8)  it  is  not  difficult  to  obtain  the  integral  rela¬ 
tions  for  the  coefficient  pm(r/,  z/)  on  the  object. 


Pm  —  /  27Tflrm(r/,  Zf,  Tq,  Zq)pm(vq,  Zq )  dil  -(-  iT 

J  n 


(14) 


where  fl  is  the  (rg,  zq)  domain  of  integration.  Once  having  solved  forpm,  the  field 
can  be  computed  at  an  arbitrary  field  point  (r/,Zf)  off  the  scattering  object  using 
the  Helmholtz  integral  of  Eq.(8);  then  the  true  scattered  pressure  value  at  this  point 
can  be  computed  by 

oo  oo 

P(»7  •*/!*)  =  E  Zi)eim6  =  E  emPm{rf,  Zf)  COS  1710  (15) 


m=— oo  m=0 

where  em  =  1  for  m  =  0  and  em  —  2  otherwise. 

The  computational  advantage  to  this  Fourier  decomposition  is  that  a  full  three- 
dimensional  problem  has  been  reduced  to  solving  a  sequence  of  much-smaller  two- 
dimensional  problems. 

For  this  program,  the  geometric  models  of  the  cylinders  have  been  hard-wired  into 
the  software.  At  program  runtime,  the  structural  configuration  (whether  flat  or 
hemispherical  ends)  and  dimensions  (overall  length  and  radius)  may  be  selecetd. 


4  Theoretical  Target  Strength 


Target  strength  (TS)  is  defined  [9, 10]  as  the  ratio,  on  a  decibel  scale,  of  the  acoustic 
intensity  (/,)  scattered  in  a  particular  direction  to  the  incident  acoustic  intensity  (/;) 
where  the  scattered  and  incident  intensities  are  determined  at  unit  distance  from  the 
acoustic  center  of  the  target.  Alternatively,  the  intensity  parameters  may  be  replaced 
by  expressions  for  pressure  (p)  under  the  far-field  approximation  which  leads  to  the 
following  expression  for  target  strength. 


TS  =  20  log  (“)  +  20  log  (“) 


(16) 


where  p„  represents  the  scattered  pressure  and  p,  represents  the  incident  acoustic 
pressure.  Note  that  the  scattered  pressure  is  measured  at  a  specific  location  in  the 
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far  field  while  the  incident  pressure  is  that  measured  at  the  target  location.  This 
equation  also  contains  the  term  to  correct  to  the  reference  measurement  distance  of 
one  unit  of  length  (r0,  usually  lm).  This  correction,  based  on  the  assumption  of 
spherical  spreading,  is  made  to  account  for  the  relative  position  of  the  measurement 
point  with  respect  to  the  acoustic  center  of  the  target,  i.e.:  where  r  represents  the 
distance  between  the  field  point  and  the  acoustic  center  of  the  target. 

For  monostatic  target  strength,  the  measurement  field  point  or  sonar  receiver  posi¬ 
tion  lies  in  the  direction  back  towards  the  source  location  and  the  scatter  is  referred 
to  as  backscatter.  For  bistatic  target  strength,  the  measurement  point  may  lie  in  any 
direction  relative  to  the  target.  Of  some  interest  is  the  scattered  response  directly 
opposite  the  target  from  the  source,  which  is  known  as  forward  scatter. 

In  this  type  of  analysis,  the  non-dimensional  frequency  ka  is  frequently  used,  where 
(for  acoustic  wavelength  X) 


and  a  is  a  representative  dimension  of  the  target  (typically  radius  for  spheres  and 
cylinders). 

The  theoretical  solution  for  the  monostatic  target  strength  of  a  finite  cylinder  of 
radius  a  and  length  L  is  given  [11]  as: 


TS  =  10  log(^p)  (18) 

which  is  valid  for 

L2 

ka  >>  1  ,  r  >  —  (19) 

A 

where  the  source  is  assumed  to  be  broadside  to  the  cylinder.  The  analyses  in  this 
report  were  done  at  frequencies  of  2.5,  5,  and  7.5  kHz,  which  correspond  to  ka 
values  of,  respectively,  2.62,  5.24,  and  7.86.  These  ka  values  then  correspond  to 
minimum  distances  from  the  cylinder  of  6.7m,  13m,  and  20m.  The  field  point 
distance  selected  was  20m.  For  the  flat-ended  cylinder,  the  above  equation  results 
in  a  target  strengths  at  lm  of  2.2  dB  at  5  kHz,  -0.8  dB  at  2.5  kHz  and  4.0  dB  at 
7.5  kHz.  This  equation  models  cylinders  with  flat  endcaps  [12]  and,  thus,  would 
not  account  for  the  differences  due  to  the  different  endcaps.  For  cylinders  with 
high  length-to-radius  ratios,  any  differences  should  be  relatively  minor.  As  the 
cylinders  used  in  this  study  were  relatively  short,  the  equation  may  be  less  accurate 
for  the  models  with  hemispherical  endcaps.  Also,  the  lowest  frequency  used  does 
not  strictly  match  the  condition  ka  »  1. 

An  approximation  to  the  forward  scattered  target  strength  of  a  cylinder  (with  the 
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source  broadside)  is  given  [13]  as: 

2  aL\ 

—)  <20> 

which  results  in  a  FTS  of  4.4  dB  at  2.5  kHz,  10.5  dB  at  5  kHz  and  14.0  dB  for  a 
frequency  of  7.5  kHz.  As  this  equation  is  based  on  the  projected  area  of  the  cylinder, 
it  should  more  accurate  for  the  cylinder  with  flat  ends. 

For  the  flat  end  cylinder,  the  monostatic  target  strength  for  the  case  with  the  source 
in  line  with  the  end  of  the  cylinder  can  be  approximated  by  that  of  a  flat  plate  which 
is  given  in  [10]  as: 

TS  =  10  log  (21) 

which  is  valid  for  the  same  range  as  the  cylinder’s  target  strength.  In  this  equation, 
A  is  the  area  of  the  plate  facing  the  source.  This  equation  results  in  target  strengths 
of  -9.7  dB  at  2.5  kHz,  -3.7  dB  at  5  kHz,  and  -0.2  dB  at  7.5  kHz. 


FTS  =  20  log 


For  the  cylinder  with  hemispherical  ends,  the  monostatic  target  strength  for  a  source 
on  the  end  of  the  cylinder  can  be  approximated  by  that  of  a  sphere  which  is  given 
in  [10]  as: 

TS  =  10  log  1"^-  +  —tan2  \  J2(ka  sira/3)j  (22) 


which  is  valid  for 


ka  >>  1  ,  r  >  a 


(23) 


where  J\  is  a  Bessel  function  of  order  1  and  f3  is  the  bistatic  angle  between  the 
direction  of  the  incident  wave  and  the  scattered  wave  (which  is  0°  for  backscatter). 
For  all  three  frequencies,  this  results  in  a  monostatic  target  strength  of  -18.1  dB  . 
If  the  formula  is  used  to  calculate  the  forward  scattered  target  strength,  it  results 
in  values  of  -9.1  dB,  -3.5  dB,  and  -0.11  dB  for,  respectively,  2.5  kHz,  5  kHz,  and 
7.5  kHz. 


5  Numerical  Models 


As  stated  above,  two  different  cylinders  were  examined  in  this  study.  Each  had 
an  overall  length  of  2.0m  and  a  radius  of  0.25m.  As  the  study  only  involved  the 
prediction  of  rigid  target  strength,  internal  structure  was  not  modelled. 

The  AVAST  models  used  are  shown  in  Figure  1.  Two  different  models  were  used 
for  each  type  of  endcap.  For  the  2.5  kHz  case,  the  flat  end  cylinder  had  a  total  of 
4710  panels,  while  for  the  higher  frequency  runs  the  model  used  had  9288  panels. 
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Figure  1:  AVAST  Cylinder  Models 
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The  hemispherical  end  cylinder  model  had  3976  panels  for  the  2.5  kHz  run  and 
8480  panels  for  the  5  kHz  and  7.5  kHz  runs.  These  models  resulted  in  21  elements 
per  wavelength  at  2.5  kHz,  15  at  5  kHz,  and  10  at  7.5  kHz.  A  rule  of  thumb  used  in 
early  model  development  was  to  insure  at  least  six  (four-  or  three-noded)  elements 
per  acoustic  wavelength.  Early  testing  showed  this  level  of  discretization  to  be 
inadequate  and  subsequent  testing  showed  ten  to  twelve  elements  per  wavelength 
to  be  the  minumum  level  of  discretization  for  reasonable  results. 

The  target  strength  was  predicted  using  field  points  at  a  radius  of  20m  in  a  plane 
coincident  with  the  cylinder  axis  and  a  source  at  a  distance  of  50m. 

For  the  axisymmetric  program,  the  models  were  obviously  much  simpler.  The 
cylinders  with  hemispherical  endcaps  were  discretized  to  an  average  level  of  ap¬ 
proximately  14  elements  per  wavelength  and  the  cylinders  with  flat  ends  to  a  level  of 
approximately  16  elements  per  wavelength.  This  discretization  was  skewed  slightly 
towards  more  refined  endcaps.  The  maximum  number  of  elements  used  was  thus 
156  at  7.5  kHz.  The  target  strength  was  calculated  with  field  points  at  a  radius  of 
10m. 

For  both  programs,  the  sound  speed  used  was  that  of  sea  water  (1500  m/s)  and  the 
target  strength  values  were  corrected  to  a  distance  of  lm.  Both  programs  assumed 
a  free- field  environment. 

6  Results 


Target  strength  predictions  for  both  computer  codes  were  made  for  both  the  bistatic 
case  and  the  complete  monostatic  case  (where  the  monostatic  target  strength  is  cal¬ 
culated  for  every  angle  of  interest  and  plotted  on  one  curve).  Thus,  for  the  bistatic 
predictions,  a  single  source  location  was  used  with  field  points  located  on  a  circle 
around  the  target,  while  for  the  second  case,  the  source  and  receiver  were  moved 
together  around  the  target. 


6.1  Bistatic  Target  Strength 

Some  representative  bistatic  target  strength  patterns  calculated  by  AVAST  and  the 
ASBIEM  software  are  shown  in  Figures  2  to  10.  The  figures  show  the  calculated 
bistatic  target  strength  patterns  for  three  source  positions  at  the  three  different  fre¬ 
quencies  for  the  two  models.  The  scale  is  in  dB  at  lm. 

The  first  set  had  the  source  broadside  to  the  cylinders  (90°),  the  second  source  was 
20°  off  the  cylinders’  axes,  and  the  third  set  shows  the  target  strength  with  the  source 
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end-on  (90°)  to  the  cylinders.  In  all  the  polar  plots,  the  source  location  is  marked 
with  an  asterisk  (*)  and  any  theoretical  results  are  marked  with  a  diamond  (o). 

In  general,  the  figures  show  good  agreement  between  the  AVAST  results,  the  AS- 
BIEM  results,  and  the  theoretical  results  (where  available).  As  expected,  the  the¬ 
oretical  results  gave  a  better  match  for  the  cylinders  with  flat  endcaps.  Between 
the  numerical  models,  the  cylinder  with  hemispherical  endcaps  seem  to  show  better 
agreement  than  those  from  the  flat  end  cylinder.  As  well,  over  the  range  of  frequen¬ 
cies  used,  the  flat-end  cylinder  target  strength  typically  showed  more  lobes  (was 
“spikier”)  than  that  of  the  cylinder  with  hemispherical  ends.  Both  of  these  results 
may  be  due  to  the  sharp  geometric  change  at  the  flat  end  which  should  cause  sharper 
changes  in  the  target  strength  and  which  may  cause  some  problems  for  the  numeri¬ 
cal  codes  which,  in  general,  deal  better  with  slightly  curving  geometries.  Typically, 
the  agreement  did  not  degrade  with  frequency,  indicating  that  the  grid  refinement 
appeared  to  be  sufficient  even  at  the  highest  frequency. 

The  dominant  feature  of  the  target  strength  with  the  source  broadside  was  forward 
scattering  with  levels  typically  5-8  dB  higher  than  the  backscatter.  Also,  there  were 
typically  more  lobes  in  the  forward  scatter  direction.  With  the  source  at  the  20° 
angle,  the  largest  difference  was  the  lack  of  an  obvious  specular  reflection  from  the 
endcap  for  the  cylinder  with  hemispherical  ends.  At  all  three  frequencies,  the  for¬ 
ward  scatter  and  the  specular  reflection  from  the  sides  of  the  cylinder  gave  similar 
patterns  between  the  two  models,  however,  particularly  at  the  two  higher  frequen¬ 
cies,  there  was  a  large  specular  reflection  from  the  flat  endcap  which  was  not  present 
for  the  other  model.  The  worst  agreement  of  this  set  of  results  appeared  to  occur 
for  the  7.5  kHz  case  with  the  source  at  0°.  The  AVAST  results  showed  strong  side 
lobes  at  ±20°  which  did  not  appear  in  the  ASBIEM  results.  Note  that  results  for 
the  cylinder  with  hemispherical  ends  were  indistinguishable  from  those  for  a  sphere 
with  the  same  radius. 
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FLAT  ENDS  HEMI.  ENDS 

Figure  2:  Bistatic  TS  at  2.5  kHz  -  Source  Broadside  ( —  AVAST,  ■  •  •  ASBIEM) 


FLAT  ENDS 


HEMI.  ENDS 


Figure  3:  Bistatic  TS  at  5  kHz  -  Source  Broadside  ( —  AVAST,  •  ■  ■  ASBIEM) 
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FLAT°ENDS  HEMI°  ENDS 


Figure  4:  Bistatic  TS  at  7.5  kHz  -  Source  Broadside  ( —  AVAST,  ■  •  •  ASBIEM) 


FLAT  ENDS  HEMI.  ENDS 


Figure  5:  Bistatic  TS  at  2.5  kHz  -  Source  Oblique  ( —  AVAST,  •  •  •  ASBIEM) 
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FLAT  ENDS  HEMI.  ENDS 

Figure  6:  Bistatic  TS  at  5  kHz  -  Source  Oblique  ( —  AVAST,  •  ■  ■  ASBIEM) 


FLAT  ENDS  HEMI.  ENDS 

Figure  7:  Bistatic  TS  at  7.5  kHz  -  Source  Oblique  ( —  AVAST,  •  •  •  ASBIEM) 
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FLAT  ENDS  HEMI.  ENDS 


Figure  8:  Bistatic  TS  at  2.5  kHz  -  Source  on  End  ( —  AVAST,  •  •  •  ASBIEM) 


FLAT  ENDS  HEMI  ENDS 


Figure  9:  Bistatic  TS  at  5  kHz  -  Source  on  End  ( —  AVAST,  ■  •  •  ASBIEM) 
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FLAT  ENDS  HEMI.  ENDS 


Figure  10:  Bistatic  TS  at  7.5  kHz  -  Source  on  End  ( —  AVAST,  •  •  •  ASBIEM) 
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6.2  Monostatic  Target  Strength 

For  each  of  the  models  and  each  of  the  three  frequencies,  complete  monostatic 
target  strength  diagrams  were  also  calculated  and  are  shown  in  Figures  11  to  16. 
In  this  case  they  are  plotted  as  X-Y  plots  rather  than  polar  plots  with  90°  being 
broadside  and  0°  and  180°  being  endfire.  As  the  systems  were  symmetric,  only  one 
half  the  target  strength  plots  are  shown.  As  above,  theoretical  results  are  marked 
with  a  diamond  (o). 

As  can  be  seen  from  Figures  1 1  through  14,  the  AVAST  and  ASBIEM  codes  showed 
very  good  agreement  for  both  the  2.5  kHz  and  5  kHz  cases  for  both  cylinder  mod¬ 
els.  At  the  lowest  frequency,  AVAST  seemed  to  show  a  slightly  better  agreement 
with  the  theoretical  results.  At  the  highest  frequency  (Figures  15  and  16),  however, 
AVAST,  while  correctly  matching  the  broadside  and  endfire  results  and  the  general 
level,  appears  to  contain  some  anomalous  large  side  lobes  (at  ±15°  and  25°  from 
broadside)  for  the  flat  end  cylinder  and  an  oddly  deep  null  (at  ±48°  from  broadside) 
for  the  cylinder  with  hemispherical  ends.  In  past  analyses,  these  side  lobes  disap¬ 
pear  when  a  more  refined  grid  is  used.  It  is  not  clear  at  this  time  what  causes  these 
anomalous  lobes  and  it  is  being  pursued.  Note  that  the  results  appear  to  match  the 
ASBIEM  results  at  all  other  locations. 

As  in  the  previous  section,  the  results  from  the  cylinder  with  hemispherical  ends 
showed  less  variation  (from  peak  to  valley)  than  the  flat  end  cylinder.  The  flat 
end  cylinder  plots  show  a  larger  number  of  peaks  and  the  nulls  were  typically  5 
to  15  dB  lower  than  those  of  the  cylinder  with  the  hemispherical  ends.  For  the 
exact  broadside  case,  the  numerical  codes  matched  the  theoretical  results  better  for 
the  flat  end  cylinder  indicating  that  the  theoretical  values  are  best  applied  to  this 
type  of  cylinder.  Based  on  a  series  of  runs  performed  with  the  ASBIEM  code  and 
other  AVAST  results,  the  flat-end  and  hemispherical-end  cylinder  results  seem  to 
converge  as  the  cylinder  aspect  ratio  increases. 

7  Conclusions 


Two  separately  developed  numerical  codes  for  predicting  the  acoustic  target  strength 
of  underwater  structures  were  compared  using  models  of  mine-sized  cylinders  at  a 
variety  of  frequencies.  While  both  programs  use  the  boundary  integral  equation 
to  solve  the  target  strength  problem,  AVAST  uses  full  three-dimensional  models 
while  ASBIEM  deals  with  only  axisymmetric  models.  While  AVAST  is  capable  of 
solving  much  more  complex  structures,  by  its  very  nature  it  will  take  much  longer 
to  solve  those  problems  which  can  be  modelled  axisymmetrically  (using  the  same 
element  dimensions). 
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Figure  11:  Monostatic  TS  of  Flat  Cylinder  -  2.5  kHz 


Figure  12:  Monostatic  TS  of  Hemi  Cylinder  -  2.5  kHz 
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Figure  13:  Monostatic  TS  of  Flat  Cylinder  -  5  kHz 


Figure  14:  Monostatic  TS  of  Hemi  Cylinder  -  5  kHz 


18 


DREATM  2001-018 


P516761.PDF  [Page:  32  of  37] 


Figure  15:  Monostatic  TS  of  Flat  Cylinder  -  7.5  kHz 


Figure  16:  Monostatic  TS  of  Hemi  Cylinder  -  7.5  kHz 
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The  models  analysed  were  those  of  cylinders  of  overall  length  2.0m  and  radius 
0.25m.  The  models  were  examined  at  three  frequencies  (2.5  kHz,  5.0  kHz,  and 
7.5  kHz)  in  a  free-field  environment.  Results  were  generated  for  both  the  bistatic 
and  full  monostatic  cases  and  compared  to  theoretical  results  where  available. 

For  the  bistatic  case,  there  was  generally  good  agreement  between  AVAST,  AS- 
BIEM,  and  theoretical  results  for  the  specific  source  locations  selected.  The  results 
from  the  cylinder  with  hemispherical  endcaps  seem  to  show  better  agreement  than 
those  from  the  flat  end  cylinder,  likely  due  to  the  sharp  geometric  changes  with  the 
flat  end  cylinder.  Typically  all  three  frequencies  showed  good  agreement  indicat¬ 
ing  that  the  level  of  refinement  was  sufficient.  The  dominant  feature  of  the  bistatic 
target  strength  in  general  was  forward  scattering  which  was  always  the  largest  am¬ 
plitude  lobe. 

For  the  monostatic  analyses,  the  AVAST  and  ASBIEM  codes  showed  very  good 
agreement  for  both  the  2.5  kHz  and  5  kHz  cases  for  both  cylinder  models.  At  the 
lowest  frequency,  AVAST  seemed  to  show  a  slightly  better  agreement  with  the  the¬ 
oretical  results.  At  the  highest  frequency,  however,  AVAST,  while  correctly  match¬ 
ing  the  broadside  and  endfire  results  and  the  general  level,  appears  to  contain  some 
anomalous  side  lobes  and  nulls.  In  past  analyses,  these  side  lobes  disappear  when 
a  more  refined  grid  is  used,  so  while  the  bistatic  results  implied  the  grid  was  suffi¬ 
ciently  refined,  these  results  indicate  this  was  not  the  case.  At  least  for  the  AVAST 
code,  ten  elements  per  acoustic  wavelength  is  not  quite  sufficient  for  the  complete 
analysis. 

In  general,  the  above  comparisons  show  that  AVAST  and  ASBIEM  give  comparable 
results  when  used  to  predict  rigid  target  strength.  The  ASBIEM  code  runs  on  the 
order  of  one  magnitude  faster  than  the  AVAST  code,  but  can  only  be  used  to  model 
structures  which  are  axisymmetric.  Future  work  on  AVAST  will  concentrate  on 
examining  elastic  target  strength  (allowing  for  internal  structure),  exploring  ways 
to  increase  the  frequency  range  (including  examining  the  anomalous  lobes  demon¬ 
strated  here),  and  examining  a  variety  of  environmental  conditions  such  as  shallow 
water  and  varying  bottom  conditions. 
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